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BaBar Results on CP Violation in B Decays 

Romulus Godang (On Behalf of the BaBar Collaboration) 
Department of Physics, University of South Alabama, Mobile, Alabama 36688 

We report on the study of the decay B+ — > D°(D0)K + where D° and DO decaying to Ktttt , with the Atwood 
Dunictz and Soni (ADS) method. We measure the ratios Rads, R~ since the processes B+ — > D°K+ and 
B+ — > D°K+ are proportional to V c j, and V u t, respectively, are sensitive to tb and to the weak phase 7. 
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Ph 1. Introduction 

m . 

During recent years, several methods have been proposed to obtain the information on the Cabibbo- 
Kobayashi-Maskawa (CKM) d 0, phase angle 7. In the Standard Model (SM), the angle 7 is the relative 
phase between b — > cus and b — > ucs transitions as indicated in Fig. Q] and Fig. [2j 
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Figure 1: Feynman diagrams for £> + — > D°K + (b — > ucs)) 
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Figure 2: Feynman diagrams for B + — > D if + (6 — > cus) 



This phase angle 7 can be measured using a variety of methods involving P-meson decays that mediated by 
either only tree-level or both tree-level and loop-level amplitudes. A theoretical source of information on the 
angle 7 is provided by B~ — > D^K~ decays. The represents an admixture of and 5"° states. 

These decays exploit the interference between B~ — > D^°K~ and B~ — > D^°K~ that occurs when the D^*^° 
and Z)(*)° decay to common final states. 

In the Gronau-London-Wyler (GLW) method [1 [|, the D meson is reconstructed based on Cabibbo- 
suppressed to CP-eigenstates, such as K + K~. In order to determine the angle 7 from B ± decays, we define 
the two experimental observables direct-CP-violating partial decay rate asymmetries as 

_ T(B- -» D C p±K-) - T{B+ -» D C p±K+) ^ ±2r B sin6 B sin 1 
CP± ~ T{B~ -> D C p±K-) + T(B+ D C p±K+) ~ 1 + r% ± 2r B cosS B cosj ^ ' 

We also define the two ratios of charged averaged partial rates using D meson decays to CP and flavor eigenstates 
as 



(2) 
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where Dcp± refer to the CP eigenstates of D meson system and 8b is the difference of their strong phases, and 
rs is the magnitude of the ratio of the amplitudes for each decay 

_ \A(B~ ^ D°K-)\ 
TB ~ \A(B- -► D°K-)\ [ ' 

In the Adwood-Dunietz-Soni (ADS) method [!, ||, the D° meson is reconstructed in the doubly Cabibbo- 
suppressed decay D° — > K + tt~ from the favored b — > c, while the D° meson from the favored b — > u suppressed 
amplitude is reconstructed in the favored decay D° —> K + ir~ . By ignoring the possible effect due to the D 
meson mixing, we define the the charge ratio of B + and B~ decay rate to the ADS final states R + and R~ , 
respectively. 

R+ = - [ ^7 + -\Z2 =rl+r 2 D + 2r B r D k D COS ( 7 + S B + S D ) (4) 



and 



T{B+ -» [K+tt-]K+) 



T{B- [K+ir-]K-) ,2 
R = r/ r>_ r „_ +1 = r B + r D + 2rBrDKD cos(7 — Ob + on) (5) 



T(B~ [K-ir+]K 

where rs and rr> are the suppressed to favored B and D amplitude ratios as 



\A{B+ -> D°K+))\ y ' 



and 

2 r(L>° K+tx- 



and ^ b and are the strong phase differences between the two B and the two D, respectively. 



(7) 



2. BaBar Detector 

The results presented in this paper are based on the entire BB data sample collected with the BABAR detector 
at the PEP-II asymmetric-energy B factory at the SLAC National Accelerator Laboratory. The BB pairs are 
produced from the decays of T(4S) resonance (on-resonance) that originate in collisions of 9.0 GeV electrons 
and 3.1 GeV positrons. The on-resonance data sample has a mean energy of 10.58 GeV and an energy rms 
spread of 4.6 MeV. The off-resonance (continuum) data sample has a center-of-mass (CM) energy 40 MeV 
below the resonance. 

A detailed description of the BABAR detector and the algorithms used for track reconstruction and particle 
identification is provided elsewhere Q- A brief summary is given here. High- momentum particles are recon- 
structed by matching hits in the silicon vertex tracker (SVT) with track elements in the drift chamber (DCH). 
Lower momentum tracks, which do not leave signals on many wires in the DCH due to the bending induced by a 
magnetic field, are reconstructed in the SVT alone. Electrons are identified by the ratio of the track momentum 
to the associated energy deposited in the calorimeter (EMC), the transverse profile of the shower, the energy 
loss in the drift chamber, and information from a Cherenkov detector (DIRC). The BABAR detector Monte Carlo 
simulation is based on GEANT4 @|. We use EVTGEN [t| to model the kinematics of B meson decays and use 
JETSET [l0| to model off-resonance process e + e~ — > qq (q = u, d, s, or c quark). 



3. Data Sample 

In the GLW method [ll|, we use (467±5) x 10 6 BB pairs, approximately equally divided into B°B° and B + B~ . 
The data have been collected in the years from 1999 until early 2008. We have reconstructed B^ — > DK^ decays, 
with D mesons decaying to non-CP (Kir), and CP-even (K+K +, tt+tt-), and CP-odd (K°n° , K°<p, K° s uj) 
eigenstates. 

In the ADS method, we use two results where D° mesons are reconstructed into two modes, _D° — > K + n~ [l2[ 
and D° -> K + 7r~ir° [l|, respectively. In D° K + n~ mode we use a data of (467 ± 5) x 10 6 BB pairs. We 
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present a search of the decays B~ — » D^*'K~ , where the neutral D mesons decay into K + tt~ final state (WS). 
In this paper, we first applied to B~ — > D^tt" , where the neutral D mesons decay into the Cabibbo-favored 
(K~ir~) and doubly suppressed mode. In D° — > K + ir~ir°, we use (474 ± 5) x 10 6 BB pairs. An additional 
off-resonance data sample of 45 fb~ x , collected at a center-of-mass energy 40 MeV below the T(4S) resonance, 
is used to study the e + e~ — > qq background. In this paper, we studied the decays D° and D° in which decays 
to K^w^ir final states. 



4. Measurement of the CKM angle 7 

4.1. CP Observables in B ± DcpK^ (GLW) 

We identify signal B — > DK and B — > Dir candidates using two defined kinematic variables. The first 
variable is the difference between the CM energy of the B meson (Eg) and the beam energy (AE) and the 
beam-energy-substituted mass (mEs), respectively. 

AE = E* B - v/i/2 (8) 

We estimate the irreducible background yields in our sample by exploiting the fact that the D invariant mass 
distribution fir this background is approximately uniform, while the signal event peaks around the nominal D 
mass. 

Figure [3] shows the AE projections of the final fits to the CP subsamples. The curves are the full PDF (solid, 
blue) and B — > Dn (dash-dotted, green) stacked on the remaining backgrounds (dotted, purple). The region 
between the solid and the dash-dotted lines represents the contribution of B — > DK. Figure |4] shows tties 
projections as well as projections to the fit to the D° — > K~ir + flavor mode. The line definitions are the same 
as described in Fig. [3] We obtain the most precise measurements of the GLW parameters Acp± and Rcp±- 

A C p+ = 0.25 ± 0.06(stai) ± 0.02(sys<) A CP - = -0.09 ± 0.07{stat) ± 0.02(syst) (9) 



R C p+ = 1.18 ± 0.09(stat) ± 0.05(syst) R CP - = 1.07 ± 0.08(stat) ± O.OA(syst) (10) 

We measure a value of Aqp+ which is 3.6 standard deviations from zero, which constitutes the first evidence 
for direct CP violation in B —> DK decays. 




AE (GeV) AE (GeV) 



Figure 3: AE projections of the fits to the data, split into subset of define CP of the D candidate and charge of the B 
candidate: a) B~ — > Dcp+K~ and b) B + — > Dcp+K + 
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Figure 4: m^s projections of the fits to the data, split into subset of define CP of the D candidate and charge of the B 
candidate: a) B~ — ¥ Dcp+K~ and b) B + — ¥ Dcp+K + 



4.2. CP Observables in B~ DK~ and B~ D*K~ (ADS) 

We reconstruct £?~ — > D^*'K~ and £>~ — > D^ir' with the D meson decaying to if~7r + (RS = right-sign) 
and fsT + 7r~ (WS = wrong-sign). Charged conjugate reactions are assumed throughout this paper. For decays 
involving a D* , both D* — > Dir° and D* — > D7 are reconstructed. To study the BB background for each signal 
category and charge combination (RS and WS) we use a sample of e + e~ — > T(4S) — > BB Monte Carlo events 
corresponding to about 3 times the data luminosity. 

In this paper the off-resonance background events are reduced by using a multilayer perceptron artificial 
neural network with 2 hidden layers, available in the framework of TMVA package [14J. We use the neural 
network to select the discriminating variables. 

Figure [5] shows the projections on tues and neural network (NN) of the fit results for DK + mode for samples 
enriched in signal with the requirements NN > 0.94 for tues projections or 5.2725 < ities < 5.2875 GeV/c 2 . 
The point with error bars are data. The curves represent the fit projections for signal plus background (solid), 
the sum of all background components (dashed), and qq background only (dotted). The results of fits to the 
B + and B~ sample 

R+ = (2.2 ±0.9 ±0.3) x 1(T 2 FT = (0.2 ± 0.6 ± 0.2) x 10~ 2 (11) 

and we extracted the variables r B 

r B = (9.5±H)% r* B = (9.6±^)% (12) 



4.3. CP Observables in B ± [K^n ± 7z°]K ± (ADS) 

In this paper the off-resonance background events, in contrast to BB events, are characterized by a jet-like 
topology. We use a Fisher discriminant J- to discriminate between the two categories of events. The Fisher 
discriminant is a linear combination of six variables. We choose the coefficients of the linear combination 
to maximize the separation between the signal and the off-resonance background. For the signal the Fisher 
discriminant T peaks at 1 and -1 for the off-resonance events. Since the correlations among the variables are 
negligible, we write the PDFs as products of one dimensional distributions of the vcies an d T. We use MC 
samples to check the absence of the correlation between these distributions. 

The PDF parameters are derived from data sample when possible in order to reduce the systematic uncer- 
tainties. The parameters for the continuum events are determined from the off-resonance data sample. We 
use the data sample of B + — > Dir + with D — > K + Tt~-K a to extract the parameters for the wies distribution. 
The parameters for the non-peaking BB distributions and the signal Fisher PDF are determined from the MC 
sample. 
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Figure 6: Constraints on from combined B — > [Ktv]K ADS measurements (left) and constraints on angle 7 from 
combined B~ D^*' [K + ir~]K~ ADS measurements (right). 



The fits to the ities for T > 0.5 and the Fisher discriminant T distribution with rriES > 5.27 GeV/c 2 for 
the combined B + and B~ samples are shown in Fig. [7] and Fig. [SI respectively. The data arc well described 
by the overall fit results (solid blue line) which is the sum of the signal, continuum background, non-peaking 
BB background, and peaking BB background. We have presented a study of the decays B — > D°K ± and 
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Figure 7: Distribution of mss (a,b) with T > 0.5. 
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B ± D ^, in which the D° and D° mesons decay to K ± tt ± tt final state using the ADS method. The final 
results are 



R + = (5±g±i) x 10" 



R- = (12i-±i) x 10- 



(13) 
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Figure 8: Distribution of T (c,d) with uies > 5.27 GeV/c 2 . 



from which we obtain 90% probability limits 

R+ < 23 x 1CT 3 R~ < 29 x 1CT 3 (14) 

Following a Bayesian approach [ID, [lj| , the probability distributions for the R + and R~~ ratio obtained in the fit 
are translated into a probability distribution for Tb- Figure [9] shows the posterior probability distribution. Since 



x10' 




Figure 9: Bayesian posterior probability density function for rs from our measurement of R and R and the hadronic 
D decay parameters ro, Sd, and ko- 

the measurements are not statistically significant, we integrate over the positive portion of that distribution 
and obtain the upper limit re < 0.13 at 90% probability, and the range 

r B e [0.01,0.11] (15) 

at 68% probability and 0.078 as the most probable value. 
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